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Abstract 

We perform a QCD sum rule analysis for the light scalar meson a (/o(600)) with 
a tetraquark current related to the instanton picture for QCD vacuum. We demon- 
strate that instanton current, including equal weights of scalar and pseudoscalar diquark- 
antidiquarks, leads to a strong cancelation between the contributions of high dimension 
operators in the operator product expansion (OPE). Furthermore, in the case of this cur- 
rent direct instanton contributions do not spoil the sum rules. Our calculation, obtained 
from the OPE up to dimension 10 operators, gives the mass of a-meson around 780MeV. 
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1 Introduction 



Nowadays there is a lot of controversy in the interpretation of the scalar mesons with 
the masses below 1 GeV [1.. In the constituent quark models, they are expected to have 
the quark content of qq as the normal members of flavor SU (3) / nonet with one unit of 
orbital excitation for positive parity. However, from the fact that the orbital excitation 
gains energy about 0.5 GeV, it is difficult to explain their light masses as well as their 
mass spectrum (so called inverted mass spectrum) |2j. Additionally, the two candidates 
to the members of nonet, the isovector ao(980) and isoscalar / (980), have a very peculiar 
properties. Indeed, their masses are degenerated and they have strong coupling to KK 
channel in strong contradiction with expectation of simple qq picture of the mesons. 
This puzzle stimulated alternative interpretations of these mesons as various types of 
tetraquark states e.g., meson-meson molecule states 4], diquark-antidiquark bound 
states [2 13 El El El , and as some hybrid states of mixture of qq and meson-meson and 
of diquark-antidiquark and meson-meson [TU] . 

In this paper, we will consider the properties of the lightest scalar meson state a 
(/o(600)) as a diquark-antidiquark bound state within the QCD sum rule ^1] with the 
special choice of the tetraquark current. We should mention that the QCD sum rules 
(SR) for the light scalar mesons were already considered separately with the interpolating 
currents of the scalar diquark-antidiquark and of the pseudoscalar diquark-antidiquark. 
In these calculations only contributions to OPE up to the operators of dimension 6 have 
been considered [3IH]. However, it was recently shown in [TJJ [T31 [TJ| that SR for multi- 
quark systems might receive a large contributions from the operators of higher dimensions 
which can lead to strong instability of the obtained results for some types of interpolat- 
ing currents. In particularly, it has been demonstrated that in the case of the scalar 
diquark-scalar antidiquark interpolating current for light tetraquarks, the contribution 
of the operators of dimensions 8 violates the requirement of positivity of left-hand side 
(LHS) and leads to disappearance of the bound state tetraquark signal [T5] . 

In this Letter we suggest to use special type of interpolating tetraquark current for 
lightest scalar meson a which leads simultaneously to the cancelation of high dimensional 
condensate contributions to the OPE and some dangerous instanton contribution to SR. 
The cx-state has the vacuum quantum numbers and should couple to the QCD vacuum 
very strongly. Instantons, topological fluctuations of gluon fields, play very important role 
in structure of QCD vacuum [13] and in spectroscopy of the multiquark hadrons [TH], [T7j . 
|18j . Therefore, our basic idea is to use the color-spin-flavor structure of the four-quark 
interaction induced by instantons JH] to fix the possible "good" interpolating current 
for er-state. This specific interaction gives strong correlations between scalar diquark of 
3 C and scalar antidiquark of 3 C , between pseudoscalar diquark of 3 C and pseudoscalar 
antidiquark of 3 C , and between tensor diquark of 6 C and tensor antidiquark of 6 C , where 
the subscript c means color. With this inspection, we propose that the interpolating 
current for cx(600) consists of the above three types of diquark-antidiquark combinations. 
Our additional argument in favor of instanton current is based on the phenomena of the 
cancelation of high dimension operator contributions expected for the case of the OPE in 
the self-dual vacuum fields (20] • Instanton field is a self-dual field and this self-duality 
property manifests directly in the color-spin-flavor structure of the four-quarks instanton 
interaction. Therefore, we can expect a similar cancelation for the case of tetraquark 
current obtained from the quark-quark instanton induced interaction. 
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We construct the QCD SR with the OPE up to operators of dimension 10 and show 
that for the case of instanton current with equal weights for the scalar diquark-antidiquark 
and the pseudoscalar diquark-antidiquark, the cancelation takes place separately for the 
high dimension operator contributions and for some dangerous instanton contributions 
to the SR. Our results for the left-hand side (LHS) of the SR are very stable. For the 
phenomenological, right-hand side (RHS) of the SR, we apply the two resonance approx- 
imation for the spectral representation, which allows to avoid the well known problem of 
strong dependence of the results for multiquark systems on the value of threshold [T^ I21j. 

This paper is organized as follows. In Sec. II, using the instanton induced quark- 
quark interaction, we fix the general structure of interpolating current for cr-tetraquark. 
In Sec. Ill we construct the standard QCD sum rule for a based on the OPE. In Sec. IV 
the direct instanton effects to the sum rule are considered. We present the numerical 
results in Sec. V and discuss them in the Conclusion. 



2 The instanton current for cr-tetraquark 

The famous 't Hooft instanton induced interaction between light quarks [22! for case of 
Nf = 2 can be written in following form ^J] 



C 



G 



4(JV* 



2N r . 



2N 



(1) 



where ip is two flavors spinor, N c is the number of colors and r~ = (r, i). This Lagrangian 
can be transformed to the Lagrangian for the interactions between diquark and antidi- 
quark by a Fierz transformations in the spin, flavor and color spaces : 



C 



G 



. v ( y , ) (^ T C 7 V 2 A^(^ 2 AV^ T ) + (Vt2\ A C^t 2 \ a C^ 
(/ {^ T r 2 X s CcrP^)^T 2 X s a p ^ T ) 



16N C (N C + 1) 



(2) 



where \ A ' S are the antisymmetric and symmetric color generators normalized as Tr(A a A b ) = 
25 ab , respectively. The first two terms correspond to the scalar and the pseudoscalar di- 
quarks in antisymmetric 3 C representation and the last one to the tensor diquark in 
symmetric 6 C color state. 

By introducing the spin matrices, 



Ts = C^ 5 , T PS = c, r Tj 
we can rewrite the Lagrangian in terms of flavors 



pa G <J pa , 



(3) 
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where Tj = 7°r|7° and the properties 

r 5 = -c 7 5 , i 



PS 



C, Y 



T,pcr 



"CpcrC, 



(4) 



(5) 
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have been used. One can see that only a very restricted set of diquarks can strongly 
interact with instanton field, namely scalar, pseudoscalar and tensor diquarks. 
Therefore we suggest the following interpolating current for the a-meson : 

r = ar s + pj° s + 7 r T (6) 

where each current is defined by 

Js = £abceade(u[T s d c )(u d T S (%) , 

J PS — e abc^ade{u^T PS d c ){u d T ps3^) , 

J£ = (ulT TtPa d a >)(u a T^dJ, + Ua'T^d^), (7) 

where a, (5 and 7 are some constants. From the Lagrangian Eq. (fl|). for N c = 3, it is 
expected that the ratio of the coefficients 

a : $ : 7 = 1 : -1 : \ (8) 

may provide some specific properties of OPE expansion for the a correlator and may 
finally lead to the most stable QCD sum rule. 



3 The OPE contribution to cr-meson correlator 

The a-correlator for the case of current Eq. (JHJ) is decomposed into nine parts 
IF = if d 4 xe iq - x (0\TJ a (x).r j (0)\0) 

= a 2 n^ + /5 2 n P5 ' P5 + 7 2 n T ' T 

+ap(u s ' PS + u PS ' s ) + « 7 (n 5 ' T + n T < 5 ) + /3 7 (n P5 > T + n T > P5 ) . (9) 

H A ' B means the correlator between A-type current and S-type current. Full set of the 
diagrams for the cx-correlator is presented in Fig. ^ It is evident that we should consider 
only the last diagram contribution to SR since only that diagram is relevant in the de- 
scription of a as the tetraquark state. The propagator for massless quarks q = u, d in 
the fixed point gauge in Fig. ^ up to order of the g 2 in strong coupling constant is the 
following [23j 

S-Jx) = -i(0\Tg a (x)q„(0)\0) 

where a, b are the color indices. 

The OPE, up to the operators of dimension 10, yields the imaginary part of the cx- 
correlator 



-ImlWg 2 ) = (a 2 + (3 2 + A8^ (V) ' 
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(a) 



Figure 1: Diagrammatic representation of the various correlators Eq (jHJ). Each quark line 
means the full quark propagator. 
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where each term corresponds to each diagram shown in Fig |2] and the factorization hy- 
pothesis for high dimension operators has been used. 




(e-l) (e-2) (e-3) (f-1) ^ (f-2) 

Figure 2: Diagrammatic representation of Eq. (fTTj) . 



4 The direct instanton contribution 

In addition to contributions of power type from the OPE expansion to the QCD SR, there 
are exponential contributions coming from direct instantons contributions as shown in 
Fig. El Their contributions can be calculated by using the following formula in Euclidean 
space for the quark propagator on the instanton background in the regular gauge 

S q J nst (x,y) = A^yh^il + ^iUr-r+U^U, (12) 

where 

P 2 

Aq(x, y) = -i J )(x - z )<p(y - z ) 
lb7T z m* 

and 

Here p stands for the instanton size, zq for the center of the instanton. U represents the 
color orientation matrix of the instanton in SU(3) C and t^'~ are SU(2) C matrices. The 
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effective mass of quark on the instanton vacuum is m* = m q — 2n 2 p 2 (qq) /3 with current 
quark mass m q . At the final stage, we multiply the result by a factor of two to take into 
account the anti-instanton contribution and integrate over the color orientation and the 
instanton size. 




Figure 3: Direct instanton contributions to the correlator. 



With the definition Q 2 = —q 2 , the direct instanton contribution to the cr-correlator 
from the two diagrams in Fig. El is given by 



w I+ AQ) = (« 2 -/3 2 )^#/ 6 (Q) 



^effPiiqq)' 



+ [19(a 2 + (3 2 ) - Qa(3 + 912^ + 72a 7 + 72(3 7 ] \ J ^ C ™' f (Q), (13) 



187r 4 m* 



where Shuryak's instanton liquid model for QCD vacuum with density n(p) = n e ff8(p — 
pc) [IH] has been used and I means the contribution from anti-instanton. fe(Q), fo(Q) 
are the functions defined by 



HQ) = J d A z 1 d A x 



x 6 [z 2 + plfKx - z ) 2 + p 2 } 



3 



e 



tq-x 



HQ) = Jd%Jd*x [4 + pimx _ ZQ)2 + pl]3 . (14) 

Here p c is the average instanton size. Moreover, let us also note that direct instanton 
contribution is possible only for the different quark flavors. Therefore, in case of tetraquark 
a meson with pure u- and d- quark content, there is no direct instanton contribution in 
any subsystem of three quarks. On the other hand, three-body instanton contribution 
in uds,uds,uds subsystems might be quite important in the case of / (980) and a (980) 
tetraquarks. Furthermore, such interaction can also lead to sizeable mixing of udud o~- 
meson with convenient ss two-quark state. 



5 Numerical analysis of QCD sum rule for instanton 
current 

In order to avoid strong dependence of the multiquark mass on the value of threshold 21J, 
we apply the two resonances approximation to the spectral representation of the correlator 

as 

Imn ff (s 2 ) = 7iY,S(s 2 -m 2 n )(0\J a \n)(n\J^\0) 

n 

= 27if 2 m s 1 5(s 2 - ml) + 27if 2 m 8 2 S(s 2 - m\) + 6(s 2 - sg)Im U OPE \s 2 )(15) 
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with the convention 

(0| r\S { ) = V2f t mt . (16) 
In this case the QCD sum rule is the following 

2 , a2 , „ Q 2 ,M W E 4 | / 2 i ol i qq 2 i in ^ 2 G 2 >M 6 £ 2 



(a 2 + /3 2 + 48 7 2 )^J + (a 2 + ^ + 88 7 2 + 12a 7 - 12/3 7 )- ^ 
+ (« 2 - ^) T ^M 4 E 1 - (a 2 - ^) 1^ M^ 



■ 2 - (« 2 - ■ ^ 

12tt 2 V P ; 12tt 2 

,,2 2 ^9((z<7ga • Gg)) 2 2 q2 J(9 2 G 2 )(m} 2 , ,.2 2 , 32n e// p^ 

+ (« -Z 9 ) 2 1 " ■ 3 2 7T 2 + ( ~ ^ > 26 ■ 33tt2 + ^ ~^ ) ^8^2 B [M<3)] 

+ [19(a 2 + P 2 ) - 6aP + 912 7 2 + 72a 7 + 72Pj) ^0^ B[f (Q)) 
= 2flm\e~ m "l M2 + 2/ 2 2 m«e- m '/ A/2 (17) 

up to operators of dimension 10, where B[f 0;6 (Q)] is the Borel transform fofi(Q) function 
given by Eq. (jl4j) . The contribution from the continuum is encoded in the functions E n (M) 
defined by 

E n (M) = — 1 ———[ 4 ds 2 e~ s2 / M \s 2 ) n , (18) 

v ; T(n + l)M 2n + 2 Jo v ' ' K J 

where so is the threshold of the continuum and M is the Borel mass. In the calculation 
of the Borel transformed fe(Q) we only include the contribution from the pole at finite 
distance x 2 ~ — p 2 , which corresponds to the direct instantons effect (see discussion in 

B2|): 

7T 4 M 12 r 1 r 1 e -M 2 pl/{ity{l~y)) , 

BIMQ)} = -—j o Mf o Ay y2(i _ y)2 (X> + 5X* + WX> 
+ 10X 5 + 5X 6 + X 7 

B[f (Q)} = ^ e -^/ 2 (ir (MV 2 /2)+^(M^/2)) , (19) 

where X = (1 — t)/t and K n (x) is the McDonald function. 

For the numerical analysis with the massless quarks q — u,d, we use the following 
condensates at normalization point p = lGeV and the average size of instanton, 

(qq) = -(0.23 GeV) 3 , (igqa ■ Gq) = 0.8 GeV 2 (gg), 

( g 2 G 2 ) = 0.5 GeV 4 , p c = lAGeV- 1 . (20) 

As we already mentioned, the sum rules constructed with only the scalar diquark-antidiquark, 
(5 = 7 = 0, and with only pseudoscalar diquark-antidiquark, a = 7 = 0, are not sta- 
ble. More precisely, the sum rule with the scalar diquark-antidiquark looses its physical 
meaning because the LHS of the sum rule has definite negative value On the other 
hand, the LHS of the sum rule with the pseudoscalar diquark-antidiquark has a definite 
positive value but its slope is negative so that it is impossible to fit a physical mass for the 
resonance. We present this situation in Figs. 4 and 5 with the value s — 1.0 GeV which 
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Figure 4: The left hand side of the QCD sum 
rule for a with the scalar diquark and the 
scalar antidiquark and so = l.OGeV. 
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Figure 5: The left hand side of the QCD 
sum rule for o with the pseudoscalar di- 
quark and the pseudoscalar antidiquark and 
s = l.OGeV. 



is the threshold taken usually in the single resonance approximation for the RHS of SR 
EI] ■ The origin of this behavior of the LHS of the sum rule lies in large contributions 
from the higher dimensional operators and direct instantons. 

Therefore, the pure scalar or pseudoscalar diquark content of tr-tetraquark is not fa- 
vored in the QCD sum rule approach and we need to find another interpolating current 
which will be not so much affected by the higher dimension condensates and the direct 
instantons. Let us discuss in detail the origin of the specific dependence, of the different 
terms in the LHS of the SR Eq.([17|h on the parameters of current a and (3. We would 
like to emphasize, that the chirality structure of the current plays an important role in 
the appearance or disappearance of some OPE and instanton contributions 1 . Indeed, 
the inspection of our general current Eq.© shows that the two quarks (antiquarks) in 
the diquarks (antidiquarks) in all above currents should have the same chirality. The 
appearence of an extra 75 in the scalar diquark-antidiquark current induces an addi- 
tional negative sign between opposite chirality diquarks with respect to the pseudoscalar 
diquark-antidiquark current so that 

aJg+(3Jp S ~ -(a - f3)(ulCd L u L Cd[ + u R Cd R u R Cd^) 

+(a + p)(u R 1 Cd R u L CdJ + u T L Cd L u R C¥ R ) , (21) 

where we have dropped the color indices for simplicity. From Eq. (jUJ) it is evident that there 
are several possible contributions to the SR with definite chirality flip. Firstly, the chirality 
conserved diagonal transitions between four terms in Eq. (j21|) gives the factor (a 2 + f3 2 ) 
in front of chirality conserved OPE terms in the SR. Secondly the contribution is coming 
from a non-diagonal transition between the first and second line in Eq. © . In this case a 
flip of chirality of two quarks in the four quark system happens and the factor (a 2 — j3 2 ) in 
front of the high dimension OPE condensates, and some instanton contributions, appears. 
There is also the possibility to have a chirality flip for all quarks coming from transitions 

1 We should mention, that from our point of view, the chirality structure of the multiquark currents 
might be one of the cornerstones of multiquark spectroscopy because it provides a strong restriction on 
the possible "good" multiquark currents. Particularly, in recent papers ^] [13] it has been shown that 
the chirality structure of the pentaquark current is very important to explain its small width. 
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between first (third) and second (fourth) terms. In OPE expansion this contribution is 
not vanishing only for massive quarks. Indeed, one cannot put all quark lines to zero 
virtuality in chirality odd condensates due to necessity to have non-zero momentum 
transfer through the last diagram in Fig^ This momentum transfer with chirality flip 
in the corresponding quark line can be produced in the OPE only for a mass dependent 
term in the perturbative part of the full quark propagator. Therefore, for a-tetraquark 
such contribution is not possible. In addition to the chirality arguments above, the spin 
structure of the tensor current 

r T ~ -{u R Ca pu d R + u L Ca pa d L ){u R a^Cd R + u L a^Cd L ) (22) 

restricts further the high dimension OPE contributions from the tensor current to the 
sum rule. It is evident from our discussion above that the cases a = —[3 and a — (3 are 
very peculiar cases for the a SR. With these values of a and (3 the contributions from 
operators with dimension higher than d = 6 and some of direct the instanton contributions 
in Eq. (jl7j) vanish. Therefore, for these currents we could expect more stable QCD SR 
results. The difference between the two cases is again in the chirality structure of the 
corresponding current. For a = —f3 each term in the first line in Eq.([21|) has four units 
of chirality, while for a = f3 each term in second line in Eq. (j21|) has zero chirality. From 
the point of view of the instanton model for the QCD vacuum the two cases presented in 
Fig© 




a) b) 



Figure 6: Two suitable currents for a- meson a) single instanton current and b) cur- 
rent induced by instanton-antiinstanton molecules. The symbol / (/) denote instanton 
(antiinstanton). 

The FiglHK corresponds to our single instanton current which we discussed in Sec. II. 
This current could be related to the phase of instanton liquid with spontaneous chiral 
symmetry breaking [15] . The FigEb corresponds to the chirality symmetric component 
of the QCD vacuum. Within the instanton model one can consider it as contribution 
coming from instanton-antiinstanton molecules. The latter component is expected to 
give a small contribution at zero temperature but might be important at temperatures 
above the deconfinement temperature. Therefore, we will choose a = —(3 as a "good" 
interpolating current for cr-tetraquark based on the instanton picture of the QCD vacuum 
and the chirality arguments given above 2 . Assigning a = —f3 = 1,7 = 1/4 to the total 

2 We should point out, that if we will choose a = j3, the analysis of SR gives for the value of the 
tetraquark mass approximately the same as for the case of "good" current for the threshold s = 2 GeV. 
The physical meaning of the small splitting of the tetraquark masses with two completely different currents 
is difficult to explain for the present. 
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current based on the structure of the instanton Lagrangian Eq.fjHJ), we have stability in 
the sum rule thanks to the cancelation in the contributions from the higher dimensional 
operators and some part of the instanton effects. With the best fit in the two resonance 
approximation, the fitted masses and the residues are summarized in Tab. 1 and the good 
quality of the fit is shown in Fig. [7| We interpret the lower mass in our SR as the mass of 
cr-tetraquark and the higher mass as the mass of its first radial excitation. According to 
the PDG the candidates for the excited states are / (1370), / (1500), and / (1710). 
We can see that the lower mass mi almost does not depend on the value of the threshold 
and the higher mass m 2 looks similar to /o(1710). Indeed, since the quality of the fit of 
the masses and the residues is best at sq = 2.0 GeV, we can interpret the lower mass 
~ 780 MeV as the mass of the a and the higher mass ~ 1775 MeV as the mass of the 
state /o(1710) 3 . In this scheme the /o(1370) might be treated as a conventional qq state 
and the / (1500)-meson can be consider as a glueball candidate 0. 



so (GeV) 


m t (GeV) 


m 2 (GeV) 


A(10- 3 GeV) 


/ 2 (10- 3 GeV) 


2.0 


0.7822 


1.7756 


4.399 


0.395 


2.2 


0.7964 


1.9488 


4.241 


0.426 


2.4 


0.8102 


2.1016 


4.095 


0.459 


2.6 


0.8261 


2.2519 


3.935 


0.491 



Table 1: Fitted masses and residues in the two resonances approximation with a = 1, /3 = 
-1,7 = 1/4. 
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Figure 7: The left hand side (solid line) and the right hand side (dashed line) of the QCD 
sum rule with masses and residues presented in Tab. 1 at Sq — 2.0 GeV. 



6 Conclusion 

We have discussed a novel interpolating current for scalar meson <r(600) treated as a 
tetraquark state. The color-spin-flavor structure of our current is fixed by the the proper- 
ties of the instanton induced quark-quark interaction and reflects the topological structure 

3 Recently, this meson was reported to have a mass /o(1790) by BES collaboration [25] . 
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of the QCD vacuum. Our "good" current has very peculiar chirality structure because 
it includes both type of diquarks, scalar and pseudoscalar, with equal weights. In this 
connection we would like to point out that the similar improvement of the OPE conver- 
gence for some specific currents was previously observed in the case of the usual hadrons. 
Well known example is so-called Ioffe's current for the nucleon which is widely using for 
description of the nucleon properties [27]. This current has also very peculiar chirality 
structure and includes scalar and pseudoscalar diquarks with equal weights. 

We have demonstrated that the definite chirality structure of our "good" current for 
a leads to the cancelation of the high dimensional operator contributions to OPE and 
to the vanishing of some instanton contributions which can spoil the QCD sum rules. 
As a result, within QCD sum rule approach, we have obtained a very stable result for 
the mass of the a meson around 780 MeV. This mass lies inside the range of the PDG 
m <r(600) — 400-^-1200 MeV. We should also mention the possible change in the result for the 
mass if the mixing between the usual two quark states and our tetraquark state is taken 
into account. An additional shift of the a-mass might come from the large perturbative 
QCD corrections in multiquark hadrons [2*%] . 
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